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Abstract. We discuss precision measurements of the leading atmospheric parameters at a standard neutrino factory. The 
oscillation of the muon and electron neutrinos (anti-neutrinos) to tau neutrinos (anti-neutrinos) adds to the muon events 
sample (both right sign and wrong sign) via leptonic decays of the taus produced through charge-current interactions in the 
detector. We focus on how this contribution affects a precision measurement of the atmospheric mixing parameters and the 
deviation of <-> v T mixing from maximality. 
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INTRODUCTION 

A neutrino factory, with decays of ji + or beams in the 
straight sections of a storage ring, produces a spectrum of 
electron and muon neutrinos via jJ. + ~^ Vf 1 e + v e and 
V^e~V e - Charged-current (CC) interactions of the muon 
anti-neutrinos (neutrinos) in a detector lead to the pro- 
duction of ji + (p.~) or the "right-sign" (RS) muons. The 
electron neutrinos (anti-neutrinos), on the other hand, 
can oscillate into muon neutrinos (anti-neutrinos) pro- 
ducing muons with charge opposite to that of the unoscil- 
lated case and get detected as "wrong sign" (WS) muons. 

Most studies with neutrino factoriesfl] have focused 
on the use of WS events to pin down the unknown reac- 
tor angle 613, the CP violating phase 8cp and the mass 
hierarchy. RS muons are useful for precision measure- 
ments of the atmospheric mixing angle 023 an d the mass 
squared difference Am 2 , apart from understanding cross- 
section and flux uncertainties. In particular, they are sen- 
sitive probes of whether <-> V T mixing is maximal (i.e. 
#23 = 7t/4, referred to as maximal 623)- Measurement of 
deviation from maximality is of significance in develop- 
ing models for neutrino masses and mixings. 

We address the little-studied issue of contamination 
of the (RS or WS) muon events sample from oscilla- 
tions of the muon or electron neutrinos (anti-neutrinos) 
to tau neutrinos (anti-neutrinos), which, through CC in- 
teractions, result in tau leptons that decay to muons. We 
focus on how this contribution affects a precision mea- 
surement of 023 and its deviation from maximal mixing. 

The oscillation probability P^ T is large even if 613 van- 
ishes, being driven by a nearly maximal 023- In spite of 
CC cross-section suppression for the massive tau produc- 



tion, there is still a sizeable production rate and tau decay 
(with rate into muons of ~ 17%) enhances the RS muon 
event rates, especially at small muon energies. Since the 
023 dependent terms in P^ and P^ T come with opposite 
sign, the combination of muons from direct production 
and from tau decays marginally decreases the sensitivity 
of the event rates to this angle. Cuts imposed to remove 
the tau events drastically reduce the direct muon events 
as well, worsening the sensitivity to the oscillation pa- 
rameters. Thus, the total muon events (direct muons and 
from tau decay) are less sensitive to the deviation from 
maximality. Neglect of the tau contribution will lead to 
an incorrect conclusion about the precision possible for 
the deviation from maximality. 



THE INPUTS: FLUXES, KINEMATICS, 
CROSS-SECTIONS 

The neutrino factory fluxes 

We assume a basic muon storage ring configuration [2] 
with muon beam energy Ef, = 25 GeV, with n^ = 5 x 10 20 
useful decays per year. We integrate the neutrino fluxes 
over the muon beam angle a, assuming a gaussian angu- 
lar divergence of the muon beam around the z-axis with 
standard deviation [3] <7 = 0.1 //, where 7 = Eyjm^. 
This is then averaged over a small neutrino opening an- 
gle, 0' < e = 0.3<7 or roughly 0. 1 mr. The resulting neu- 
trino (or antineutrino) spectrum with gaussian angular 
spread averaged over (0' < e,0 < 0' < 27r) at the base- 
line distance L from source to detector is given by, 
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The kinematics 

We focus on the spectrum of the final state muons and 
hence require the detailed kinematics of the CC interac- 
tions, in which either muons or taus are produced, with 
the latter decaying into muons, where again we use the 
differential decay rates (see Ref. [4] for details). In the 
laboratory frame, a neutrino of flavor / = ji or x interacts 
with a nucleon and produces the corresponding charged 
lepton / at an angle 0/ w.r.t. the incident neutrino direc- 
tion. In the case of v T interactions, the tau is produced 
at a very forward angle while the azimuthal angle of the 
muon from tau decay 0^ , is restricted by the decay kine- 
matics. The available phase space is restricted in both 
direct muon and tau production due to the constraint on 
the available energy for the lepton: £_ < E\ < E + ; see 
Ref. [4] for the detailed expressions on the constraints. 
The effect of this pinching in available energy, for the 
case of a T lepton being produced, can be seen in Fig. 1 
where the final hadronic mass mx is plotted as a function 
of E T . The notation is standard: m x — W 2 — (p + q) 2 , 
where p, q are the nucleon and intermediate gauge bo- 
son 4-momenta in the laboratory frame. For a typical 
neutrino energy E v = 10 GeV, allowed energy range for 
cos T =0.91 is tiny; tau leptons are hence produced in 
a very forward direction while the direct muons, due to 
their lighter mass, are less restricted. 



The cross-sections 

Since the energies of interest range from a few GeV to 
25 GeV, the CC interactions include quasi-elastic (QE), 
resonance (Res) and deep inelastic (DIS) processes. We 
consider the double differential cross-sections, 
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where Gp is the Fermi constant, k — Mfy / (Q 2 is 
the propagator factor with the W boson mass, Mw, Pi is 
the magnitude of three-momentum of the charged lepton 
produced and W, are structure functions corresponding 



FIGURE 1. Kinematics of v T -nucleon CC interactions. The 
allowed parabolas of constant cos0 T in the mx-E? plane are 
shown for E v = 10 GeV. The ends of the parabolas (at mx = M) 
give the limits of the tau energy. 



to the general decomposition of the hadronic tensor. W4 5 
appear only for massive final leptons. We have, 

£ a/Wi = \2W } + (E, - pi cos 6,) + W 2 {E, + p, cos &,) 

± -± (e v E, +pf- (E v +E,)p, cos 6,) ~W S . (3) 

The detailed expressions for W, are taken from Ref. [5] 
where the specific structure functions are listed for QE, 
Res and DIS leading order (LO) processes. 



EVENT RATES IN A FAR-DETECTOR 
Preliminaries 

We assume that the neutrinos interact with a 50 kton 
iron detector such as the proposed INO/ICAL or MIND. 
Both and beams with equal exposure are con- 
sidered. While RS events are sensitive to deviations of 
023 from maximality, inclusion of WS events may only 
marginally worsen the results; however, the advantage 
in being "charge-blind" is significant, hence, all muon 
events are simply added. The generic number of muon 
events at a distance L, as a function of the observed muon 
energy E is, 
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FIGURE 2. Muon event rates as a function of the observed 
muon energy. RS and WS events from fi~ and fl + beams are 
shown in the upper panels. Contributions from direct muon 
production (denoted by D) and that of muons from tau decay 
(labeled as t) are separately shown. The left lower panel shows 
the sum D + z. Oscillation parameters are as given in text. 



where the superscript b refers to the beam-type, i = 
jj.,e(jl,e) denotes the neutrino type from beam. 
Here K — N t n y , N t is number of target nucleons (as- 
sumed isoscalar), n y is number of years of data and 
R is the Gaussian energy resolution function of width 
<T = O.OlEn. The subscripts D and T correspond to the 
production of direct muons and T production and decay 
into muons, respectively. The limits of integration and 
the restriction from an angular constraint, in case of tau 
decay are defined in Ref. [4]. The charge-blind events 
are obtained by summing over the index i and the to- 
tal events for each beam are then obtained by the sum, 
0% ot {E) = Zi(^i lD (E)+^(E)) . We assume 90% re- 
construction efficiency of muons. Muon event rates accu- 
mulated over five years are used to study the sensitivity 
to the deviation of the mixing angle 023 from maximality. 

Typical event rates at L=7400 Km (magic baseline) for 
oscillation parameters, Am 2 = 2.4 x 10~ 3 eV 2 , 023 = 42°, 
013 = 1°, sin 2 12 = 0.304 and A 2 i = 7.65 x 10~ 5 eV 2 , 
are shown as a function of the observed lepton energy in 
Fig. 2. The panels show the direct muon production and 
tau decay contributions to the RS and WS muon events 
from and beams. Since the neutrino-nucleon 
cross-sections are larger than those for anti-neutrinos, the 
(/i + ) beams are preferred for studying the RS (WS) 
events where the cross-sections are larger. However, the 
full sample without charge identification is preferable, 
due to higher detection efficiency. We therefore add the 
events (RS+WS) from both and ju+ beams. 

It can be seen that there is a substantial contribution 



FIGURE 3. Effects of angular cuts on the tau contribution to 
muon events at neutrino factories. For details, see the text. 



to the RS events from tau decay into muons. The tau 
contribution alters the sensitivity to the oscillation pa- 
rameters; we next discuss whether these events can be 
removed through suitable cuts. 

Cuts on tau contribution 

Since tau production in neutrino-nucleon interactions 
is extremely forward-peaked, one obvious way to remove 
the tau contribution is with an angular cut. Also, since the 
tau contribution is substantial at small observed muon 
energies where the tau decay rate is large, a muon energy 
cut can also be contemplated. The effect of cuts on the 
event rates is seen in Fig. 3 - the only cut effective in 
removing the tau contribution is one (0^ > 25°) that re- 
moves the signal itself! Alternately a muon energy cut of 
E > 10 — 15 GeV can substantially remove the tau con- 
tribution, still leaving sufficient direct muons. However, 
such a large energy cut will worsen the measured preci- 
sion of the mixing parameters as sensitivity is higher in 
the lower energy bins where matter effects are large. In 
short, it is not feasible to cut out the tau contribution and 
still make a precision measurement, in this case, of the 
deviation of 023 from maximality. 



Effect of the tau contribution 

As stated earlier, the tau contribution alters the depen- 
dence on the mixing parameters, altering the precision 
to which we can determine them. While the tau events 
have less sensitivity to 023, the rate increases while the 
direct event rate decreases, as 023 increases towards max- 
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imal 023 = tc/A. The inclusion of muons from tau events 
also alters the uncertainties considerably. A near detec- 
tor sensitive to muons, measures the combination of flux 
times cross-section of the muons. This also appears in the 
RS event rate for direct muon production and is there- 
fore well constrained. However, for tau production and 
decay, the RS event rate has the combination of muon 
flux and the tau production cross-section. The heavy tau 
cross-sections have larger uncertainties, where mass cor- 
rections are large. Furthermore, since these contributions 
result from oscillations, no near detector can help reduce 
the uncertainties. Hence overall uncertainties are much 
larger for the tau contribution than for direct muons. 

Hence, in our numerical calculations we use an over- 
all normalization error of 0. 1 % for direct, while a mod- 
est 2% is used for the total (direct+tau), muon events. 
We use typical input values of (Am 2 , 023, 613) to estimate 
how well the generated "data" can be fitted, and calculate 
the resulting precision on the parameters. We keep the 
solar parameters fixed at their best-fit values of Ref. [6] 
and set 8cp to zero. The best fits (and regions of confi- 
dence levels in parameter space) are obtained by mini- 
mizing the chi-squared with a pull corresponding to the 
normalization uncertainties specified. 

We present results for a typical sample set of input 
parameters, (Am 2 , 23 , 9 13 = 2.4 x 1(T 3 eV 2 , 41.9°, 1°). 
We minimize first over the pull, and then over Am 2 and 
023, keeping 613 fixed. Fig. 4 shows the allowed Am 2 - 
023 parameter space at 99% CL. The solid and dashed 
lines correspond to considering direct and total (includ- 
ing those from tau decay) muon events respectively. Note 
that the 99% CL contour is much more constrained with 
direct than for total muons. In particular, it is the Am 2 
values that are smaller than the input value, that broaden 
the contour and limit the discrimination. The largest true 
value of 023 that can be discriminated from maximal is 
shown in Fig. 5, as a function of Am 2 again, for 0i3 = 1°. 
It is seen that tau contamination worsens the ability to 
discriminate 623 from maximal, thus making this mea- 
surement harder than originally expected. 



CONCLUSION 

The oscillations of the muon or electron neutrinos (anti- 
neutrinos) to tau neutrinos (anti-neutrinos) results in tau 
leptons produced through CC interactions in the detector 
which on leptonic decay add to the right as well as wrong 
sign muon events obtained directly. This tau contamina- 
tion worsens the ability to discriminate against maximal 
Vn <-> V T mixing. It is practically impossible to devise 
satisfactory cuts to remove this tau contamination. Un- 
certainties from this tau background must be brought un- 
der control before making precision parameter measure- 
ments at neutrino factories. 
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FIGURE 4. Allowed Am 2 -6 2 3 parameter space at 99% CL 
from CC muon events, directly produced (solid line) and with 
inclusion of those from tau decay (dashed line). See text for 
values of input parameters. 



2.6 



2.0 



Exclusion plot at 99%CL: Fixed 8 13 =1 





43 



44 45 46 
(Input) 9^3 



47 48 49 



FIGURE 5. The largest (smallest) true value of 2 3 in the 
first (second) octant that can be discriminated from 2 3 = ft/A, 
as a function of Am 2 are shown when only Direct (D) and total 
(D + t) events are considered; Here 6*13 is fixed at 1°. 
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